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Abstract 

We present in detail a theoretical model for fragmentation reactions of three-body 
halo nuclei. The different reaction mechanisms corresponding to the different pro- 
cesses are described and discussed. Coulomb and nuclear interactions are simul- 
taneously included and the method is therefore applicable for any target, light, 
intermediate and heavy. Absolute values of many differential cross sections are then 
available as function of beam energy and target. We apply the method to frag- 
mentation of 6 He and n Li on C, Cu and Pb. A large variety of observables, cross 
sections and momentum distributions, are computed. In almost all cases we obtain 
good agreement with the available experimental data. 

PACS: 25.60.-t, 25.10,+s, 21.45. +v 



1 Introduction 



The peculiar structure of halo nuclei is clearly revealed by fragmentation 
reactions [1-4]. Large interaction cross sections and narrow momentum dis- 
tributions of the fragments observed after fragmentation of these nuclei are 
direct signs of their large spatial extension. The main properties of these nu- 
clei are well described by few-body models [5,6], where the nucleus is viewed 
as an inert core surrounded by a few (usually one or two) bound nucleons. 
Among them, Borromean two-neutron halo nuclei have attracted a lot of at- 
tention, and their most prominent examples, 6 He and n Li, have been widely 
investigated. 
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During the last decade a large amount of experimental data after fragmen- 
tation of 6 He and 11 Li on different targets have been provided. In particular 
total dissociation cross sections and total interaction cross sections, as well 
as momentum distribution of the fragments after dissociation and after core 
breakup are available [7-18]. Different models describing these data have been 
developed. These models fall in two independent main groups according to the 
kind of reaction for which they were designed. The first one considers nuclear 
breakup processes, and therefore they apply to the case of light targets. Most 
of them describe neutron dissociation reactions while core breakup is usually 
omitted [19,20]. Only calculations performed in the framework of the Glauber 
theory are also providing total interaction cross sections [21-23]. The second 
group focuses on heavy targets where only Coulomb dissociation is considered 
[24-27]. Intermediate targets where nuclear and Coulomb interaction compete 
is outside the scope of these models. At best Coulomb and nuclear breakup 
are computed in independent models and subsequently simply added [28,29]. 

The reaction mechanisms responsible for the breakup processes are still 
controversial [30-32], but absolutely essential to understand and analyze the 
data properly [33]. It is established that few-body aspects of the halos are 
necessary. In most investigations halo constituents (and targets) are treated 
as inert particles where the intrinsic degrees of freedom are neglected or ac- 
counted for by finite range potentials with Pauli forbidden states. However, 
the core degrees of freedom are unavoidable for some observables like inter- 
action cross sections where core destruction is substantial. Such many-body 
aspects must be included at some point while maintaining the necessary halo 
features. 

More than one reaction mechanism is necessary in a correct description of 
halo breakup. The long range Coulomb and the short range nuclear potentials 
already lead to clear-cut differences [30,31]. However the finite extension of 
halo particles and targets furthermore produces inevitable differences in the 
active reaction mechanisms at small and large impact parameters. This has 
strong implications for the understanding and the analyses. 

The starting point for developing our model was the recognition that, for 
large beam energies (<; 100 MeV/nucleon), the time scales allow application 
of the sudden approximation for two-neutron [34,35] and one-neutron halos 
[36]. Only light targets, relative momentum distributions and the specifically 
dominating processes were accessible within the sudden approximation. The 
next step was to use a participant-spectator separation and employ the phe- 
nomenological optical model to describe the interactions between participants 
and target [37]. Again only light targets were accessible but now absolute cross 
sections were computed. The model was then extended to treat the interac- 
tions between spectators and target in the black sphere model [20]. Finally 
the Coulomb and nuclear interactions were included on the same footing al- 
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lowing treatment of light, intermediate and heavy targets and computation of 
absolute values of all differential breakup cross sections [30,31]. 

The reaction mechanisms now come out as clearly distinguishable for light 
and heavy targets or rather for small and large impact parameters. The small 
impact parameters involve violent collisions and one or more of the halo par- 
ticles do not arrive in the detectors in the forward direction. At large impact 
parameters the momentum transfer is relatively small and the three halo par- 
ticles are gently excited into the continuum state, which subsequently decays 
and all halo particles most likely arrive in the forward detectors. These two 
extremes are separated roughly at impact parameters slightly larger than cor- 
responding to grazing collisions. These different reaction mechanisms imply 
that the two-step processes, where one halo particle is removed by the target 
while the resonance of the remaining system is populated followed by decay 
[14-18], at best only is consistent with small momentum transfer [33]. The 
large momentum transfer requires a completely different mechanism of halo 
particle removal. 

In this paper we formulate for the first time in detail the model accounting 
for all these effects. So far no other model contains all the necessary fea- 
tures. The idea is to describe the fragmentation reaction as the sum of all 
the possible reactions where one, two or all three halo constituents interact 
simultaneously with the target. The processes leading to different final states 
clearly add up incoherently. The cross sections for the remaining processes are 
in all the numerical applications in this report also simply added as appropri- 
ate for large momentum transfer from the target. The model is consistently 
describing two-neutron halo fragmentation on any target with simultaneous 
treatment of nuclear and Coulomb interactions. Absolute values are computed 
for all possible dissociation cross sections distinguished according to the par- 
ticles left in the final state. Two-neutron removal cross sections, core breakup 
cross sections, interaction cross sections and momentum distributions are com- 
puted as function of target and beam energy above 50 MeV/nucleon. We shall 
present results for a large variety of these observables arising from fragmen- 
tation reactions of 6 He and n Li on C, Cu, and Pb. These three targets have 
been chosen as typical examples for light, intermediate and heavy targets dis- 
playing different reaction mechanisms due to the interplay of Coulomb and 
nuclear interactions. Details of the model are given in Sections 2 and 3. Re- 
sults for fragmentation reactions of 6 He and n Li on C, Cu and Pb are shown 
in Section 4. Finally Section 5 contains summary and conclusions. 
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Fig. 1. Sketch of the reaction and coordinates used for the final state two by two 
separation. The target is labeled by and {i,j,k} label the particles within the 
three-body projectile. 

2 Model and method for inert particles 

We consider a spatially extended three-body halo nucleus colliding with 
a target at high energy. Let us first assume that the target and the three 
halo constituents in the projectile are inert particles. The differential cross 
section da is to a good approximation given by the sum of three terms da^\ 
each of them describing the independent contribution to the process from the 
interaction between the target and the halo particle %. This is the assumption 
used in the classical formulation for a weakly bound projectile [38]. We neglect 
the binding energy of the initial three-body bound state compared to the high 
energy of the beam. 

The reaction is then described as three particles independently interacting 
with the target as if each particle were free. Each individual interaction with 
the target is viewed as removal of one particle % (participant) while the other 
two particles j and k (spectators) both survive the reaction undisturbed (see 
Fig. 1). The participant is either absorbed (in the optical model sense) or 
elastically scattered (diffracted) by the target. 

The masses, coordinates and conjugate momenta are denoted m, r and p, 
respectively. The three halo particles and the target are labeled by {i, j, k} = 
{1,2,3} and 0, respectively. The relative coordinates rjk, 7*oi, fi,jk, R and 
R' are defined as shown in Fig. 1. The corresponding conjugate momenta 
are analogously denoted by p jk , p 0i , p ijk , P = p Gijk and P' = p' 0i>jk . We 
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use primes to denote the momenta in the final state. Detailed expressions 
for the different relative coordinates and momenta can be found in [20]. The 
momentum transfer in the reaction for elastic scattering of the participant is 
then given by 

q = P -p'o = Poi-Poi » 1 = \q\ = 2poiSin- , (1) 



where 9 is the angle between p 0i and p' 0i . 



2.1 Two by two separation in the final state 



For relatively large momentum transfer between participant and target 
the final state is appropriately described as two independent subsystems, i.e. 
the target plus participant and the two spectators. This separation is not 
a matter of convenience but dictated by physics, especially in connection 
with the generalization in the next section to spatially extended particles. 
This is a central assumption reflecting the reaction mechanism. It is discussed 
in previous publications and tested by comparison with experimental results 
[12,14,16,17,33,34]. The transition amplitude of the process in Fig. 1 can then 
within the mixed plane/distorted wave Born approximation be written as 
[20,39] 

= (^^.^^l^l^e^) , (2) 

fQi i trjk jk jk 



where J and M are the total angular momentum of the projectile and its 
projection on the z-axis, s'j k and E^ fe are the spin and its third component of 
the two-body system formed by particles j and k after the collision, and E^ 
are the spin and projection quantum numbers of particle % (for convenience we 
assume a spin zero target), and an d 4>ry~) s' are the distorted wave 

t^Oi i ir jk jk jk 

P' Ft' 

functions of the two final two-body systems. We use plane waves e and 
P Ft 

e because these relative momenta are much higher and the plane wave 
approximation is suitable. 

In the frame of the projectile (p i + Pj + p k = 0), and assuming that the 
participant i has spin or 1/2, the previous transition amplitude leads to the 
following differential diffraction cross section for the process in Fig. 1 

d 9 ^(P',P' jk , g) d^Hpu -> p'J PM {JM) 

dP'dp' 3k dq dq 2J + 1 Mg ^ s , |M ^^I , W 
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MP'S E' = (4'1> V j P ^ ri >*Xs>V.\^ JM) ) 

jk jk i * jk jk jk i i 



(4) 



Note that the final two-body wave function 0^,* 2, is contained in the differen- 

tial elastic cross section for the participant-target scattering d 3 a^ jdq. The 
two-body wave function 0«/ fc ~/ E / is computed as in [35] and the appropriate 

Mr j k jk jk 

phase factor is included in the computation of the overlap matrix element. 
When the participant % in the sense of the optical model is absorbed by the 
target the corresponding differential absorption cross section factorizes as 

dP'd P ' 3k ~ {m) 27TT M /£ E , 1 M '»W 1 • (5) 

jk jk i 



The final state wave function used in eq.(2) is not orthogonal to the initial 
wave function. This means that our final state contains a non-zero component 
of three-body bound state plus target, that represents elastic scattering of 
the halo nucleus as a whole, and has to be removed. This is done in eq.(3) 
by means of the function P dis (q), that is the probability of dissociation of the 
three-body halo system. For elastic scattering the final state wave function 
of the halo nucleus is e %( i^' ri ^ k ^!^ JM \ where q cm = q(rrij + m^jivrti + rrij + 
m^) is the momentum transfer into the participant-spectators relative motion 
described by the coordinate r it j k . Therefore the probability of elastic scattering 
of the three-body projectile as a whole is given by the square of the overlap 
between the initial and final three-body wave functions. This implies that the 
probability of dissociation is given by [40] 

P dis {q) = 1 - |(^(^)]e^- r ^|*( JM ))| 2 , (6) 



which vanishes as g 2 (oc q 2 cm ) for small values of the momentum transfer. 

In eqs.(3) and (5) draf^ jdq and a a ° b l ] are the differential elastic and ab- 
sorption cross sections for the participant-target scattering, respectively [41]. 
We use a central participant-target nuclear potential Viv(r), and assuming 
that the participant and the target have charges Zi and Z , respectively, the 
participant-target interaction is 

V 0i (r) = ^ + V N (r) , (7) 



where a = e 2 /Tic is the electromagnetic coupling constant. 

The cross sections arising from the interaction V 0i in eq.(7) are obtained 
after integration of eqs.(5) and (3) when particle % is absorbed and elastically 
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Fig. 2. Sketch of the reaction and coordinates used for the final state three by one 
separation. The notation is as in Fig. 1. 

scattered, respectively. 

2.2 Three by one separation in the final state 

The scheme shown in Fig. 1 arises when the halo breakup is caused by 
a sufficiently violent and fast participant-target reaction. However this is not 
true in processes in which contributions from small momentum transfers are 
dominating. In these cases the small kick of the participant excites the halo 
rather gently and the breakup process proceeds via the created wave packet. 
An obvious example is halo breakup on a heavy target caused by the Coulomb 
participant-target interaction, which clearly is a dominating process, espe- 
cially at low beam energies [30]. 

This type of gentle breakup reaction mechanism requires that the final state 
is that component of the three-body continuum halo wave function which at 
large distances approaches plane waves with relative momenta specified by 
p' jk and p'ij k - The corresponding reaction is sketched in Fig. 2. The transition 
matrix in eq.(2) should then be replaced by 



r(0 = <<>U X-Jej Xs> A e tP '- R \V 0l (r 0t )\¥+\R, r jk , r iJk )) (8) 



where R now is the relative coordinate between the three-body system and 
the target, and P and P' are its conjugate momenta respectively in the initial 
and final state - note that in section 2.1 P' is the conjugate momentum of 
the relative radial coordinate between the center of mass of the two final two- 
body systems (see Fig. 1), while now P' is the conjugate momentum of the 
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relative radial coordinate between the final three-body system and the target. 
The two spin wave functions \ are labeled by their respective final state total 
and projection quantum numbers (s^) and (sj k H'j k ) defined earlier. 

The distorted wave \I/p^ satisfies the equation 

4 +) = 0, (9) 



Ti 2 d 2 
'2MdR 2 



+ Voiiroi) + H(r jk , r iJk ) - E 



where M is the projectile-target reduced mass and H(rj k , r it j k ) is the intrinsic 
hamiltonian of the projectile. A gentle kick by a long range Coulomb potential 
implies that the halo absorbs the transfered momentum as a whole, meaning 
that the adiabatic approximation should be adequate for this reaction type. 

We therefore as in [42] replace the intrinsic hamiltonian H(rj k , r^) with 
a constant, namely the ground state energy — B of the halo. This ground state 
energy, typically a fraction of an MeV, is much smaller than the total energy 
E, typically of the order of hundred MeV, and can therefore be neglected. 
The internal coordinates rj k , r it j k are now no longer dynamical variables in 
the equation but parameters and therefore dH = dr 0i leading to 



h 2 d 2 

'2Mdvl 



+ V 0i (r 0i ) - E 



^ +) = 0. (10) 



The solution to this equation is 

4 +) =^\r 0t )^(r jk ,r hjk ), (11) 

where V , p"^( r oi) is the distorted wave describing scattering of a particle with 
mass M in the potential V 0i and the function tp(rj k ,r it j k ) is determined from 
the large r 0i asymptotics 

+/(«)— U^ry*). (12) 
roi J 



Indeed, the plane-wave term of this expression describes a projectile in the 
ground state \I/( JM ) moving with the momentum P against the target, i.e. 



e 



lPr ^(r 3k ,v hjk ) = e lP - R ¥ JM \ (13) 



which gives 
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where \ii = (rrij + m k )/{mi + rrij + m k ). The matrix element T"W i n eq. (8) 
then becomes 



x K w (r w )|V-i, +) (r 0i )e« p "«*l J ">), (15) 
which factorizes as 



x(e iP '--|^(r 0i )|V'P +) (roO), (16) 



where q cm = /^q and q = P' — P. The second factor is the matrix element 
of the elastic scattering on the potential Voi of a particle with mass M - 
that is the projectile-target scattering with the potential Vo«. 



In addition to the breakup amplitude we obtain the elastic scattering of the 
halo as a whole simply by substituting the ground state wave function \1/( JM ) 
instead of the continuum wave function in the transition matrix element in 
eq.(8). This then leads to 

T e ( ; } = T e f } (¥ JM) \e iq ^ r ^^(JM)^ _ ( 17 ) 



For a two-neutron halo the largest contribution to the three-by-one final 
state is given at large impact parameters by the long-range Coulomb inter- 
action between the charged core and the target. Considering only this main 
contribution we can write the breakup cross-section (with % designating the 
core) as 



dP dp' jk dq dq 

* u ^ *- M,E',s' ,S' 

where the spin of the core is or 1/2 and d 3 cr^l h is simply the Rutherford 
projectile-target cross section. 

When the two-body interactions between halo fragments in the final state 
are neglected in eq.(18) the matrix element simply reduces to the Fourier 
transform of the three-body bound state wave function with the momentum 
Pijk shifted by 

Mem - 
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2.2.1 The continuum three-body wave function 

The three-body continuum wave function „, in eq.(18) has to be 

Mr jf$ ^li.jk 

computed. We first introduce the Jacobi coordinates (see e.g. [5]) x and y and 
their conjugate momenta k x and k y such that k x ■ x = r' jk ■ p' jk and k y • y = 
Pi,jk ' r i,jk- The corresponding hyperspherical coordinates p, (l are then also 
available. We omit here the label i specifying which set of Jacobi coordinates. 
We then define the continuum three-body wave function ^^kj, ( x > y) with 
the asymptotic behavior of a plane wave in six dimensions and an incoming 
hyperspherical wave 

4i (x, y) - exp(ik s • x + iky ■ y) + f*(Q K , n f M ~l Kp \ (19) 



where we have introduced the 6-momentum k = {k^, k^} and the correspond- 
ing hyperangular variables k = ^Jk% + and The total energy is then 
given by k as E = h 2 K 2 / (2m). 

We shall now express the function in terms of our adiabatic hy- 

perspherical basis functions $jMn(p, see e.g. [6,34,43]. The basis functions 
QjMn are labeled according to their asymptotic behavior, which means ap- 
proaching the corresponding hyperspheric harmonics, i.e. QjMn{p, ^) —> Yq(J1) 
for large p. Here we use the notation Yq(Q) for a hyperspheric harmonic cou- 
pled with the appropriate spin functions to the total angular momentum, 
where Q = {K,l x ,l y , L, s x , s y , S, J, M} denotes the full set of hyperradial 
quantum numbers, where s x = Sjk and s y = Sj. 

We now first shift the dependence of the function upon the five hyperangles 
fl K of the momentum k into the basis Y^(fl K ) (£ = {K, l x , l y , L}) of the ordinary 
non-spin coupled hyperspheric harmonic. Thus we get 



EiM*>y)*7(n«) 



> ( 2 °) 



where the function ^(x, y) now can be expanded in terms of our complete 
set of angular solutions ^ JM n{pM) to the Faddeev equations i.e. 

^(x, y) Xs , s , X s' A = p- 5/2 E Fi q L<J>)*JMn(p, n) , (21) 

JMn 



where q = {£, s'j k , S^, s-, £•} is the set of quantum numbers equivalent to Q 
in the decoupled spin basis. 
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Inserting expansion eq.(21) into the Schrodinger equation we obtain, after 
multiplication by {$jMn(p,ty\ an d subsequent angular integration, our usual 
hyperradial equations for Fj 9 ^ n (p) [6,34,43] 



H 3 1" Vnn - "77— ^JMn\P) 



dp 2 p' 1 L 2m 

= £ f-2P„„^ - Qnn') ^(P). (22) 



The asymptotic form eq.(19) can be expanded in terms of the hyperspheric 
harmonics Y^(Q K ) using 

exp(ik,x + ik,y) = ^^i K J K+2 ( K p)Y£(n K )Y ( (n) , (23) 



which implies the asymptotic behavior 



+/« ^^r p) )^(")n'w» ( 24 ) 

Then the asymptotic boundary condition for the hyperradial function Fj q ^ n (p) 
is found from eq.(24) 

p-WFSQniP) - E 

x(YQlnX^X^. fc >. (25) 

This boundary condition is equivalent to the scattering boundary condition 
used in [27]. 

The transition matrix element in eq.(18) reduces when q cm ■ r^jk «C 1 to 

Hf- • (V,p: , x<=j K, fc |* (JM) > , (26) 

which after squaring and summation over magnetic quantum numbers gives 
the usual factorized expression in terms of the B(E1) dipole strength function 
and the cross section [27]. Before expansion this matrix element also included 
all the higher multipole transitions. They could in principle arise from the 
nuclear as well as the Coulomb interaction. 



(2tt) 3 i K exp(-wp) 
(rep) 2 jK +^P) d & + /«' ^572 
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2.3 Momentum cutoff 



When the participant is charged the elastic cross section eq.(18) for q — > 
is dominated by the Coulomb cross section, which after integration over the 
directions of q diverges like g -3 . Since the square of the matrix element in 
eq.(18) vanishes as q 2 for small q, the integration of eq.(18) over q leads 
to a logarithmic divergence of this elastic cross section. However, to pro- 
duce dissociation, the energy transferred from target to participant (5E = 



yPo + rriQ — yp'o + m§) must be larger than the three-body binding energy 
B. When p and q = p — p' are parallel SE is maximized. For small B com- 



pared to the target rest mass 5E = B implies that qc ~ q^c = By I + rn^c 2 jp\, 

which reduces to B/v in the non-relativistic limit. Thus q must be larger than 
qL to produce dissociation. 

Another formal divergence arises for the Coulomb interaction at large im- 
pact parameters (low momentum transfer), where adiabatic motion only al- 
lows virtual excitations excluding dissociation. If the collision time At is long 
compared to the orbital period of motion T inside the projectile no transfer 
of energy occurs. These processes are excluded by the adiabatic cutoff. The 
decisive parameter is [44] 



where (5 — v/c, 7 = (1 — (3 2 )~ l l 2 and hcu a is the adiabatic cutoff energy, 
i.e. the maximum value of the equivalent photon energy for a given impact 
parameter b estimated in a sharp cutoff model [44]. Then %uj a m B, where 
B is the three-body binding energy, should be a reasonable estimate as the 
lowest energy needed for breakup. However, the softness of the electromagnetic 
excitation modes could influence this sharp cutoff value, but the effect on the 
cross sections is still only logarithmic. We leave the precise choice of fiuj a for 
adjustments in the numerical computations. 

In previous publications [30,31] we used Tiu a = B ps /n, where B ps is the 
binding energy between participant and the system consisting of the two spec- 
tators, i.e., Bp S = B — B 2s , where B 2s is the two-body binding energy of the 
two spectators. In a Borromean system B 2s is negative. This is a misleading 
connection for a slow reaction where the full projectile is adiabatically excited. 

From [45] the following classical relation between the momentum transfer 
q and the impact parameter b is derived: 





Q = 




(28) 
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where p is the momentum of the projectile and E^ in its kinetic energy. 

The adiabatic cutoff condition £ a = 1 and use of eqs.(27) and (28) leads 
then to q > q a , where 



Therefore integration of eq.(18) has to include only values of the momentum 
transfer q larger than the largest of qi and q a . The divergence then disappears. 



3 Model and method for spatially extended particles 

The two different final state divisions described in Section 2 are introduced 
for physics reasons and not for convenience. The underlying idea is that dom- 
inating reaction contributions arise from (optical model type of) absorption, 
which implies that the final state consists of the spectators and the participant- 
target system, and elastic core-target Coulomb scattering at low momentum 
transfer leading to a final state containing all three halo particles. This divi- 
sion is dictated by the reaction mechanism and seems to be unavoidable in a 
few-body treatment where the core degrees of freedom are neglected. Further- 
more, the finite extension of the projectile constituents and the target also 
requires a similar division according to the number of projectile constituents 
simultaneously colliding with the target. 

3.1 Reaction scenarios 

Let us extend the applicability of the inert particle model and incorporate 
effects of the core degrees of freedom and the finite sizes of the constituents. 
First we assume short-range interactions between each of the halo constituents 
and the target. Then we start by describing the individual constituent-target 
interaction in the black disk model, according to which the target absorbs 
the constituent inside a cylinder with the axis along the beam direction and 
otherwise leaves untouched. The radius of this cylinder is related to the range 
of the interaction and approximately equal to the sum of target and constituent 
radii. For the collision between the halo nucleus and the target we have then 
the four possibilities shown in Fig. 3: 

(a) Only one of the halo constituents, the participant, is inside its cylinder. The 
other two constituents are mere spectators, and survive untouched after the 



ZqZiC 2 hu, 



a 



(29) 




13 




Fig. 3. Scheme of the possible scenarios for the collision between the three-body 
projectile and the target. The notation indicates a halo consisting of two neutrons 
and a core. 

collision. The three constituents of the halo projectile give rise to three 
contributions of this kind. 

(b) Two of the halo constituents are participants. The third one is spectator. 
Again we have three possible reactions of this kind. 

(c) All the three halo constituents are inside their corresponding cylinders. 

(d) All the three halo constituents are outside their cylinders. For short-range 
interactions this means that the whole projectile is untouched by the target, 
and this process would only contribute to the elastic cross section. 

The scheme in Fig. 3 is only directly valid for short-range interactions and 
therefore not applicable for charged halos. This is a severe deficiency as even a 
two-neutron halo nucleus colliding with a heavy target involves a core-target 
Coulomb interaction that by far can dominate the cross section. 

To incorporate long-range interactions into the model we do the following: 
In the reactions shown in Figs. 3a, 3b, and 3c the participant-target interac- 
tion contains two parts, the full short-range interaction and the long-range 
potential related to small impact parameters (large momentum transfer). The 
large impact parameter part of the long-range potentials (small momentum 
transfer) corresponds to reactions as the one shown in Fig. 3d. 

In the particular case of two-neutron halo fragmentation this means that 
the momentum transfer between target and core is divided into large and 
small values corresponding to short and long impact parameters, respectively. 
Large momentum transfer contributions to the cross sections correspond to the 
processes shown in Figs. 3a, 3b, and 3c. They contain the nuclear core-target 
interaction and the large momentum transfer part of the Coulomb core-target 
interaction. They are described in Section 2.1, and the integral of eq.(3) over 
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the momentum transfer involves values of q larger than a certain value q g (q g < 
q < oo). Low momentum transfer contributions correspond to the process 
shown in Fig. 3d, where none of the halo constituents are inside their cylinders, 
and only the core-target Coulomb interaction contributes. Cross sections are 
obtained after integration of eq.(18), where the momentum transfer varies 
between q min = max{q L , q a } and q g . 

Low and large impact parameters are divided by the value b g = R$ + R c + 
na/2, see [44], where Rq and R c are charge root mean square radii of the 
target and the core, and a = hcZ Z c e 2 /2E kin is half the distance of closest 
core-target approach. Use of eq.(28) and the fact that p/Ekin = (7 + l)/7/3c 
then leads to the value q g of the momentum separating low (q < q g , b > b g ) 
and large momentum transfer (g > q g , b < b g ), i.e. 

* = ■ (30) 



3.2 Cross sections 



According to the number of particles surviving in the final state we have 
the possible cross sections: a , cr n , a c , a nni a nc and o" nrac , where <7 means no 
particles in the final state, and the indexes c and n refer to the core and the 
neutrons surviving in the final state. If some of the particles are absent in the 
final state they have been absorbed by the target in the optical model sense. 

The probabilities for finding the core and one halo neutron inside their 
cylinders are denoted P c and P n , respectively. These probabilities can be 
split in two parts, the probability of being absorbed by the target and the 
probability of being elastically scattered by the target: P c = P c ab + Pf and 

p pab 1 pel 

r n — r n T r n . 

When the three-body halo projectile collides with the target, each of the 
reactions shown in Fig. 3a occur with the probability Pj(l — Pj)(l — Pk), where 
the index % refers to the participant. In the same way each of the reactions 
shown in Fig. 3b and Fig. 3c occur with the probabilities — Pk) and 

PiPjPk- In this way, the total probability of a given process is obtained by 
adding the probabilities of all the processes in Fig. 3 producing precisely the 
corresponding group of halo constituents in the final state. 

For instance, the probability of a process in which only the core survives 
in the final state (processes giving rise to o~ c ) contains two terms, P^ b P^ b (l — 
P c ), that comes from one of the reactions shown in Fig. 3b in which the two 
neutrons inside their cylinders are absorbed by the target, and PfP^ b P^ b , that 
comes from the reaction shown in Fig. 3c in which the two halo neutrons are 
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Table 1 

Contributions to the different cross sections from the reactions shown in Figs. 3a, 
3b and 3c, where the nuclear interaction between the halo constituents plus the 
large momentum transfer (low impact parameter) part of the Coulomb core-target 
interactions are included. Two crosses mean that the contribution must be included 
twice due to the possible exchange of the two halo neutrons. The cross section 
&nnc contains also the contribution from the low momentum transfer (large impact 
parameter) part of the Coulomb core-target interaction sketched in Fig. 3d. 



Reaction 



Fig. 3a 



Fig.3b 



Fig. 3c 



Probability 
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absorbed and the core is scattered. 



In Table 1 we indicate which reactions in Fig. 3 contribute to the various 
cross sections. The projectile contains two neutrons and some of the reactions 
should be counted twice, since the symmetric combination also contributes the 
same amount. The probabilities shown in Table 1 are not explicitly computed. 
Instead we compute the cross sections corresponding to each of the different 
processes. For example to compute the first process in Table 1 (Fig. 3a) we 
note that Pf corresponds to the probability for elastic scattering of the core 
on the target with the related cross section given by eq.(3). The constraints 
that the two neutrons must be outside the cylinders are imposed by removing 
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appropriate parts of the initial three-body wave function, see section 4.2. All 
other processes in Fig. 3a are analogously computed by replacing the first 
probability by the related cross section and removing parts of the three-body 
wave function in the overlap in eq.(4) to account for the other probability 
factors, see section 4.2. When the probability factors appear symmetrically as 
for processes corresponding to Figs. 3b and 3c a specific choice must be made. 
The details for each of the possible cases are described in sections 3.3, 3.4 and 
3.5. 

Differential cross sections are given in eqs.(3) and (5), and total cross 
sections are obtained after integration over all the variables. The cross section 
for a process in which a certain group of constituents survive in the final state 
is then computed by summation of the partial cross sections obtained from 
each of the contributing reactions shown in Table 1. We emphasize that on 
top of the contributions in the last column of Table 1 the cross section a nnc 
also contains the low momentum transfer part due to the Coulomb interaction 
that corresponds to the process in Fig. 3d. 

From the cross sections uq, a n , a c , a nn , a nc and a nnc it is easy to obtain 
several cross sections of specific interest, as the two neutron removal cross 
section (a^n), that is given by all the processes where the core survives in 
the final state; the core removal cross section (o"- c ), given by all the processes 
in which the core is destroyed by the target; and the interaction cross section 
(07), defined as the cross section given by the reactions where the projectile 
looses at least one of its constituents: 

0"-2n = 0" c + One + &nnc , &-c = + °n + O nn , <7j = (J _ 2 n + &-c ■ (31) 

Adding the probabilities shown in Table 1 contributing to each of the 
cross sections we obtain the probabilities for the two-neutron removal and 
core breakup processes as 

P(<T-2n) = Pf + P n (l - Pc) + Pn(l - P n )(l ~ P c ) , (32) 
P(<7_ c ) = Pf . (33) 

Therefore the cross sections o"„ 2ri and cr_ c can be computed as indicated in 
eq.(31), or, alternatively, by computing the cross sections corresponding to the 
probabilities in eqs.(32) and (33). 

From eq.(32) we see that the two neutron removal cross section has three 
contributions. The first one, cross section corresponding to Pf, is the elas- 
tic scattering of the core by the target, and its differential cross section is 
given by eq.(3). Since Pf is the probability for the core being inside the cylin- 
der, only large momentum transfer, q > q g , should be included in eq.(3). 
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However, the two-neutron removal cross section also receives a contribution 
from the reaction in Fig. 3d, which contains the low momentum transfer part 
{Qmin < Q < Qg) of the core-target interaction. Thus, q is actually only re- 
stricted by q > q m in- in addition <7_ 2n also contains the contribution from the 
strong interaction between the halo neutrons and the target with the core as 
spectator. These contributions are given by the cross sections corresponding 
to the second and third terms in eq.(32). They are computed from eqs.(3) and 
(5) with the corresponding optical model neutron-target cross section with 
the overlap in eq.(4) containing the appropriate parts of the three-body wave 
function. Core breakup cross sections can be computed, according to eq.(33), 
by considering only the core-target potential and subsequent integration of 
eq.(5). 



3.3 Two-spectators contributions (Fig. 3a) 



The contributions in Fig. 3a correspond to processes in which one halo 
constituent (participant) interacts with the target, while the other two (spec- 
tators) remain untouched. They contribute to the cross sections a nn , o nc and 
a nnc as shown in the upper part of Table 1. The general form for the proba- 
bility of these reactions is Pj(l — Pj)(l — P&), where Pj is the probability for 
the halo constituent i being the participant. 

The procedure is then as discussed in Section 3.2. The participant-target 
interaction is described by the potential eq.(7) where the nuclear part is given 
by an optical potential, that takes into account both absorption and elastic 
scattering of the participant by the target. The differential cross sections are 
given by eqs.(3) and (5), and total cross sections are obtained after integration 
over all the variables. These integrations must take into account the fact that 
impact parameters for constituents j and k are larger than the radii of their 
cylinders. This is done (see section 4.2) by removing the appropriate part of 
the three-body wave function in eq.(4). 

When the core is participant the interference between Coulomb and nuclear 
core-target interactions is included in the calculation through the differential 
elastic participant-target cross section in eq.(3). It is important to keep in 
mind that the momentum transfer q in eq.(3) must be larger than q g in eq.(30). 
This restriction includes only the low impact parameter part of the Coulomb 
interaction and allows therefore precisely the processes in question. 
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3.4 One-spectator contributions (Fig. 3b) 

The contributions shown in Fig. 3b correspond to processes in which two of 
the halo constituents are inside their cylinders, while the third one is spectator. 
In principle each of the two participant constituents could be either absorbed 
or scattered by the target, and therefore the reactions shown in Fig. 3b can 
contribute to the cross sections a n , a c , a nn , a nc and a nnc as seen in the central 
part of Table 1. 

The probability for occurrence of one of the reactions in Fig. 3b is PiP 3 (1 — 
Pfc), where % and j refer to the two participants. To compute the cross section 
corresponding to this reaction we first select the constituent i or j where 
the interaction with the target is most conveniently described by the optical 
model. If possible we choose the core. Then the Coulomb interaction with 
the target and the interference between nuclear and Coulomb interactions are 
treated carefully. 

Cross sections are then computed as a process in which one of the partici- 
pants (i) interacts with the target through the corresponding optical potential. 
Differential cross sections are then given by eqs.(3) and (5) for elastic scatter- 
ing and absorption of constituent % by the target. Again total cross sections 
are computed after integration of (3) and (5) over all the variables. When 
the core is participant only the low impact parameter part of the Coulomb 
interaction is allowed, and q > q g in eq.(3). 

In this case (reactions in Fig. 3b) there is a second constituent, constituent 
j, inside its cylinder that therefore also interacts with the target. This fact is 
included in the calculations by removing the appropriate part of the three- 
body wave function in eq.(4), see sect. 4. 2. We then ensure that constituent j 
is close enough to the target to interact with it. As a result of this interaction 
constituent j can be either scattered or absorbed by the target. We then divide 
each of the cross sections computed as described above from eqs.(3) and (5) 
in two parts, one of them corresponding to absorption of constituent j by 
the target and the other one corresponding to scattering of constituent j. The 
weight of each of these parts can be obtained in the optical model for two-body 
systems, i.e. simply as a ( £j (a^ + a^) for absorption and a ( £j for 
scattering of constituent j. Here a^ s and a^ a are the absorption and elastic 
cross sections obtained by use of an optical potential describing the interaction 
between constituent j and the target. 

Finally, an additional requirement is needed when computing cross sec- 
tions. This is due to the fact that particle k is outside its cylinder, and there- 
fore not interacting with the target. This is again included in the calculation 
by removing the appropriate part of the three-body wave function in eq.(4) 
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as described in section 4.2. 

The procedure for constituents i, j and k described above permits to cal- 
culate cross sections for all the processes in Fig. 3b. It also permits to compute 
separately cross sections for processes where both, constituents % and j, are 
either absorbed or scattered, and also when one of them is scattered and the 
other one absorbed. 



3.5 Zero-spectators contribution (Fig. 3c) 

The contributions shown in Fig. 3c correspond to processes where all three 
constituents are inside the corresponding cylinders. The probability for this 
reaction is PiPjPk, and any of the three constituents can be either absorbed 
or scattered by the target. Thus, this reaction contributes to cross sections 
where zero, one, two, or three constituents survive in the final state, as shown 
in the lower part of Table 1. 

In analogy to the reactions in Fig. 3b we describe the core-target inter- 
action with an optical nuclear potential plus the Coulomb interaction eq.(7). 
We then compute cross sections from eqs.(3) and (5) as a process in which the 
interaction between projectile and target takes place through the core-target 
potential. Once more, only the low impact parameter part of the Coulomb 
interaction is allowed, and q > q g in eq.(3). 

Now, reactions in Fig. 3c, constituents j and k are also inside their cylin- 
ders, and therefore both interact with the target. As before, this is ensured 
by including in the overlap in eq.(4) the appropriate part of the three-body 
wave function (see sect. 4. 2). Again these two constituents can be either ab- 
sorbed or scattered by the target. As in the previous section, we then divide 
the computed cross sections given by eqs.(3) and (5) in three parts, corre- 
sponding to absorption of both neutrons, elastic scattering of both neutrons, 
and absorption of one of them and scattering of the other one. The weight 
of each of these three contributions is obtained in the optical model for two 
body systems, and the part of two absorbed neutrons has a weight given by 
the square of a^ s / (a^l + cr^), the part of two scattered neutrons has weight 
given by the square of cr^/(a^l + o^), and finally the part of one neutron 
absorbed and one scattered has a weight of 2a^ s a^/ (er^2 + °~eia) 2 ■ Again cr^] 
and <7g™2 are the absorption and elastic cross sections obtained by use of an 
optical potential describing the neutron-target interaction. 

This procedure permits then computations of the cross sections for the 
processes in Fig. 3c where any of the three constituents can be either absorbed 
or scattered by the target. 
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4 6 He and n Li on, C, Cu and Pb 



In this section we shall apply the method to fragmentation reactions of the 
two-neutron halo nuclei 6 He and n Li on three different targets, one heavy (Pb) 
with dominating Coulomb contribution one light (C) with dominating nuclear 
contribution and an intermediate mass (Cu) with comparable Coulomb and 
nuclear contributions. We shall first briefly give the specifications leading to 
the initial wave function. Second we describe how to choose the interactions 
between target and halo constituents. The last two subsections contain the ab- 
solute values of the computed cross sections and the momentum distributions, 
respectively. 



4-1 Three-body wave functions and intrinsic halo interactions 

The three-body halo wave functions are obtained by solving the Faddeev 
equations in coordinate space. This is done by use of the adiabatic hyperspher- 
ical approach. The Faddeev equations are written in hyperspherical coordi- 
nates, and separated into angular and radial parts. The eigenvalues obtained 
from the angular part enter into the radial equation as an effective poten- 
tial. This procedure has previously been successfully used, see for instance 
[6,34,43]. 

The neutron-neutron interaction is given in [43]. In the neutron-core in- 
teraction for 6He we include s, p and d waves. The parameters are adjusted 
to reproduce the corresponding phase shifts from zero to 20 MeV. The values 
are given in [27]. 

For n Li we include s and p-waves, and assume spin zero for the core. The 
parameters for the neutron-core interaction are adjusted to reproduce the 
available experimental data on 10 Li, in particular the presence of a low-lying 
p-resonance at 500 ± 60 keV with a width of 400 ± 60 keV [46-48], and an 
even lower lying uncertain virtual s-state at 0.15 ± 0.15 MeV [47,48]. These 
observed states cannot be used directly in computations assuming zero spin of 
the 9 Li-core. The reason is that the spin splitting, inevitably arising from the 
finite core spin of 3/2, produces two levels both of s and p-wave character. The 
measurements refer to the actual (spin-split) levels. Therefore in the present 
simplified model of zero core spin we must aim at statistically averaged s and 
p-wave energies, which then necessarily must be above the measured values. 
We have chosen the realistic model with spin splitting reproducing the avail- 
able data and simply reduced the spin splitting parameter to zero. The more 
realistic model with finite core spin is necessary for observables as the two- 
body invariant mass spectrum, but in most cases the differences are negligibly 
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small [49]. The resulting values for the parameters correspond to potential IV 
in [20]. 

It is well known that the use of realistic two-body potentials to describe a 
three-body system is leading to three-body binding energies too small com- 
pared to the experimental values. This is a general problem for few-nucleon 
systems [50]. For 6 He and n Li the underbinding is around 0.5 MeV and 150 
keV, respectively. These deviations are substantial compared to the three-body 
binding energies, but very small compared to the strengths of the two-body po- 
tentials. To recover the experimental value of the binding energy we introduce 
a phenomenological three-body interaction. The idea is that this three-body 
force should account for those polarizations of the particles, which are beyond 
that described by the two-body interactions. Thus, this three-body interaction 
must be of short range in the hyperradius p, since it only contributes when all 
three particles interact simultaneously. We choose again the gaussian radial 
shape V 3 b(r) = V 3 exp (— p 2 /6g), and the parameters for 6 He and n Li are given 
in [27] and [20], respectively. 

Inclusion of these effective three-body interactions together with the neutron- 
neutron and neutron-core interactions described above give a binding energy 
and a root mean square radius of 1.0 MeV and 2.50 fm for 6 He, and 0.3 MeV 
and 3.35 fm for n Li. These numbers are for both nuclei consistent with the 
experimental values, that are 973.4 ± 1.0 keV [51] and 2.57 ± 0.10 fm [52] for 
6 He and 295 ± 35 keV [53] and 3.1 ± 0.3 fm [54] for n Li. One has to realize 
that what is commonly called experimental radius for a halo nucleus is de- 
duced from reaction cross section measurements, and the connection between 
the few-body wave function and the cross section is model dependent [23]. 
The real test is therefore the direct comparison between computed and exper- 
imental reaction cross sections. The computed 6 He wave function corresponds 
to a probability of roughly 90% relative neutron-core p-waves and about 10% 
s-waves, while the corresponding numbers for n Li roughly are 60% relative 
neutron-core s-waves and 40% p-waves. 



4-2 Neutron-target and core-target optical potentials 

For the neutron-target interaction we use non-relativistic optical poten- 
tials obtained from relativistic potentials through a reduction of the Dirac 
equation into a Schrodinger-like equation [55]. These phenomenological poten- 
tials in the Schrodinger equation produce the same scattering data as obtained 
by use of the relativistic potentials in the Dirac equation [56]. In particular, 
we use the energy-dependent, A-independent parameterizations EDAI-C12 
and EDAI-Pb208 in [57] for carbon and lead targets, respectively. For copper 
target we use the energy and A-dependent potential EDAD-fitl also given in 



22 



Table 2 

Parameters for the core-target interaction. See eqs.(34) - (38). 



Reaction 


4 He+C 


4 He+Cu 


4 He+Pb 


9 Li+C 


9 Li+Cu 


9 Li+Pb 


r v (fm) 


1.245 


1.245 


1.245 


1.336 


1.336 


1.336 


r w (fm) 


1.570 


1.570 


1.570 


1.870 


1.870 


1.870 


a v (fm) 


0.751 


0.745 


0.733 


0.900 


0.900 


0.900 


a w (fm) 


0.651 


0.635 


0.607 


0.750 


0.750 


0.750 


a (MeV) 


70 


110 


108 


120 


220 


300 


ai (MeV) 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 




-0.008 


-0.010 


-0.014 


-0.008 


-0.010 


-0.014 


b (MeV) 


17 


21 


27 


28 


28 


29 


61 (MeV) 


-1.706 


-1.706 


-1.706 


-1.375 


-1.375 


-1.375 


62 


0.007 


0.007 


0.006 


0.0025 


0.0025 


0.006 



[57]. In the three cases the potentials fit the scattering data in the nucleon 
energy range from 20 to 1040 MeV. 

For the core-nucleus optical potential we use the parameterization given 
in [58] . The general shape of this potential is: 

U(r) = -Vf(r) - iWg{r) + V Coul (34) 

with form factors of Woods-Saxon type 

/(r) = (l+exp[(r-r^ 1 / 3 )K])- 1 , (35) 
g(r) = (1 + exp [(r - r w A^) / aJO' 1 , (36) 

and where Vc ui is the Coulomb interaction, and A the mass number of the 
target. The energy dependence of the potential is contained in the parameters 
V and W as 

V(A,Z,E) = a + ai ZA- 1/3 + a 2 E , (37) 
W(A,Z,E) = b + b 1 ZA 1/3 + b 2 E , (38) 

where Z is the charge of the target and E is the energy of the projectile. 

The parameters used for the different projectiles and targets are given 
in Table 2. For a 6 He projectile the radii, r v and r w , and the diffuseness 
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parameters, a v and (Xyj , £1X6 directly taken from [58]. The strength parameters 
in eqs.(37) and (38) are also from [58] except for which is reduced to allow 
for a large energy variation while simultaneously reproducing the experimental 
data in [22]. For a n Li projectile the parameters r v , r w , a v and a w are from 
[59], while the parameters in eqs.(37) and (38) are similar to those of 6 He but 
slightly modified to reproduce the available experimental data for 9 Li-nucleus 
scattering [10]. 

Together with the optical potential we need to specify the value of the 
distance R^ u \ which is the limiting impact parameter determining whether 
halo constituent j is participant or spectator. These cutoff parameters play 
an essential role in the cross sections. If the interaction between constituent % 
and the target is described by the optical model, the impact parameters for 
j and k in eqs. (3) and (5) must be limited in accordance with the different 
geometries in Fig. 3. This is done in practice by including in the overlap 
eq.(4) only the part of the initial three-body wave function \1/( JM ) such that 
the distances between particles j-i and k-i are either larger or smaller than 
the cutoff distances R^ u \ and R^ t depending on the reaction in question. We 
are therefore forcing particles j and k to be inside or outside a sphere around 
particle i, instead of a cylinder along the beam direction. To do the latter is 
technically much more difficult, and the use of spheres is a good approximation 
providing big advantages from the computational point of view. 

Several procedures can be followed to obtain these cutoff parameters: 

(a) Directly determined from the sizes of the constituent and the target: 

R { ! = r \[aJ 3 + A 2 J\ ifi = ^rlAT + R\ , (39) 

where the indices c and n refer to the core and the halo neutron, A t and 
A c are the mass numbers of target and core, and R n is the neutron radius. 
The value of the parameter tq may vary between 1.1 fm and 1.3 fm. 

(b) In the black disk model the absorption cross section for a particle j hitting a 
target is given by nR^ . We can then determine the cutoff radius from mea- 
sured absorption cross sections. For the projectiles and targets considered 
in this work some of these experimental data can be found in [10,22,60]. 

(c) The cutoff radii can also be determined from the measured mean square radii 
of target ((r 2 ) t ) and constituent ((r 2 )j), which can be found for instance in 
[61]. The relation to the cutoff radius is then given by 

\R C % = (r 2 )t + (r 2 ), + 2 fm 2 , (40) 

where 2 fm 2 is the square of the range of the nucleon interaction. This 
procedure gives values for R^J t very similar to the impact parameter b g 
used in eq.(30) to separate high and low impact parameters. 
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Table 3 

The cutoff parameters R^J t determined according to the procedures (a), (b) and 
(c), and the values used in the actual computations. 



Target 


Method 


RfJt (fm) 


Kut (fm) 


Rent (fm) 


C 


(a) 


2.7-3.1 


3.4-4.7 


3.7-5.2 




(b) 


2.6-3.0 


~ 4.0 


~ 4.7 




(c) 


3.9 


4.2 


4.8 




used 


3.5 


4.0 


4.8 


Cu 


(a) 


4.5-5.3 


5.2-7.2 


5.4-7.6 




(b) 


4.9-5.6 


- 6.5 


~ 8.5 




(c) 


5.5 


5.7 


6.1 




used 


5.4 


5.7 


6.1 


Pb 


(a) 


6.6-7.8 


7.4-10.4 


7.6-10.6 




(b) 


7.4-7.8 


~ 9.0 


~ 10.5 




(c) 


7.4 


7.6 


7.9 




used 


7.4 


7.6 


7.9 



These three procedures determine the range of variation for R£l t . The 
choice of one value or another is a matter of taste, and within the allowed 
variations these parameters may be used for fine tuning the cross sections. 
The different values obtained are shown in Table 3 together with the ones 
actually used in the calculations. We see that the range of variation of the 
cutoff is rather small for the neutron, and its value consequently is rather 
precisely determined. The same happens for 4 He and 9 Li on a carbon target, 
while for heavier targets like Cu, and especially for Pb, the cutoff radii vary 
more. 

The adiabatic cutoff energy in eq.(29) is chosen as huj a = 1 MeV for 6 He 
and 0.15 MeV for 11 Li, corresponding to the three-body binding energy and 
half of it, respectively. This means that n Li in the sharp cutoff model breaks 
up corresponding to virtual photon energies down to half of the binding en- 
ergy. The agreement with measured values discussed in the next section then 
indicates a connection with the soft structure revealed in the low lying dipole 
strength [17,24,27,28]. The high energy tail of the virtual photon spectrum 
must efficiently produce breakup of n Li [44]. 
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4-3 Cross sections 

In Fig. 4 we show two-neutron removal cross sections obtained after frag- 
mentation of 6 He and n Li on C (left part), Cu (central part) and Pb (right 
part). In the three cases the long-dashed line shows the contribution to the 
cross section from elastic scattering of the core by the target (core partici- 
pant), while the short-dashed line is the contribution from processes where 
the core is spectator. 

For a carbon target the core participant contribution is only significant at 
low beam energies, where it contributes by up to 25% of the total cross section 
for the 6 He projectile, and somewhat less for n Li. For large beam energies this 
contribution is at most 10% of the total. The main contribution is given by 
processes where the core is a spectator, and only the halo neutrons interact 
with the target. This contribution is shown by the short-dashed line in the 
figure as computed from the second and third terms in eq.(32). These two 
terms contain contributions from the reactions in Fig. 3a, where one of the 
neutrons and the core are spectators, and from the reaction in Fig. 3b, where 
only the core is a spectator. For comparison we also show in the figure the 
two-neutron removal cross section obtained only from the reactions in Fig. 3a 
where one neutron and the core are spectators (dot-dashed curve). 

In the central part we show the results for a Cu target. Due to the Coulomb 
interaction the core participant contribution dominates at low beam ener- 
gies, but decreases below that of the core spectator with increasing energy. 
For the 6 He projectile this crossing appears at a beam energy of around 600 
MeV/nucleon clearly larger than the 100 MeV/nucleon for n Li. Therefore, for 
energies from 100 MeV/nucleon to 600 MeV/nucleon the main contribution 
to °~-2n comes from core participant reactions for 6 He projectile and from core 
spectator reactions for n Li projectile. In other words, a Cu target resembles 
a light target for the n Li projectile and a heavy target for 6 He. This is due to 
the larger size of n Li compared to the one of 6 He. 

The general behavior of o"_ 2n for heavy targets is seen in the right hand 
side of Fig. 4, where the results for a Pb target are shown. Due to the large 
Coulomb interaction the core participant contribution is clearly dominating, 
see the dot-dashed line corresponding to the reaction in Fig. 3d, where only 
the large impact parameter part (q < q g ) of the Coulomb interaction is in- 
cluded. This contribution is for both projectiles very similar to the one coming 
from core participant (dashed), that includes the full Coulomb core-target in- 
teraction, the nuclear core-target interaction, and the interference between 
them. The contribution given by the neutron-target interaction (core specta- 
tor, dotted) is negligible at low beam energies for both projectiles, while at 
large beam energies it contributes around 15% of the total for 6 He and 35% 
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Fig. 4. Total two-neutron removal cross section, eq.(31), after fragmentation of 6 He 
and n Li on carbon (left), copper (middle) and lead (right) as function of beam 
energy. The long-dashed and short-dashed lines are the contributions to the to- 
tal from participant core and core spectator, respectively. In the carbon case the 
dot-dashed line is the contribution assuming that the core and one of the halo 
neutrons are spectators. In the lead case the dot-dashed line is the large impact 
parameter contribution of the Coulomb interaction (Fig. 3d). For 11 Li on Pb the 
thin solid line is the total two-neutron removal cross section for fiuj a = 0.3 MeV (see 
text). The experimental data are from [7,8,10,14,17]. 



of the total for 11 Li. The agreement with the measured values is overall rather 
good considering the incompatibility of some of the data points. As mentioned 
at the end of the previous section, the adiabatic cutoff hu> a is chosen to be 0.15 
MeV for the n Li projectile. The result of the expected logarithmic dependence 
on huj a is shown in Fig. 4 for the two-neutron removal cross section for the 
Pb target where the effect is largest. The thin line shows the computed cross 
section with Tiu a = 0.3 MeV. The agreement with the experimental points is 
only slightly worse. For the C target the minimum value of the momentum 
transfer is given by qi instead of q a (see sect. 2.3), and for Cu target qi be- 
comes larger than q a for beam energies around 150 MeV/nucleon when hu a = 
0.15 MeV and around 400 MeV when hu a = 0.3 MeV. Therefore the choice of 
one value or the other for the adiabatic cutoff energy does not play any role 
for C target and insignificant for the Cu target. For 6 He projectile we used 
hu a = 1 MeV, corresponding to its binding energy. A variation of this number 
within a reasonable range produces a smaller change than in the case of n Li 
due to its smaller charge. 
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Fig. 5. Total core breakup cross sections, eq.(31), after fragmentation of 6 He and 
11 Li on carbon (left), copper (middle) and lead (right) as function of beam energy. 
The long-dashed, dot-dashed and short-dashed lines are the contributions to the 
total from the reactions in Fig. 3a (two neutrons spectators), Fig. 3b (one neutron 
spectator) and Fig. 3c (no spectators). 

In Fig. 5 we show core breakup cross sections after fragmentation of 6 He 
and n Li on C (left), Cu (middle) and Pb (right). For the three targets we 
have plotted the contributions from the reactions in Fig. 3a, where the core 
is absorbed while the two neutrons are spectators (long-dashed line), the 
contributions from the reactions in Fig. 3b where only one of the halo neutrons 
is spectator (dot-dashed curve), and the contribution from the reaction in 
Fig. 3c, where there are no spectators (short-dashed line). 

For a carbon target we see that for the n Li projectile the contribution 
from the two-neutron spectators is dominant, and the contribution from si- 
multaneous collisions with the target of more than one constituents is smaller, 
although still significant. For the 6 He projectile the situation is different, and 
the no spectators contribution is even larger than the two-neutron spectators 
contribution. This is a reflection of the difference in size between these pro- 
jectiles. If we denote the neutron-core distance r nc we obtain (r 2 ^) 1 / 2 = 4.2 
fm for 6 He and (r^) 1 / 2 = 5.9 fm for n Li. In 6 He the two halo neutrons are 
then closer to the core than in n Li. Thus, when the core is absorbed the two 
halo neutrons have a larger probability of being spectators for n Li than for 
the 6 He projectile. 

For a copper target the no spectators contribution is dominant for core 
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Fig. 6. Total interaction (solid), two-neutron removal (long-dashed) and core 
breakup (short-dashed) cross sections after fragmentation of 6 He and n Li on carbon 
(left), copper (middle) and lead (right). For n Li on Pb the thin solid and dotted 
lines are the interaction and two-neutron removal cross section for huj a = 0.3 MeV 
(see text). Experimental data are from [7,8,10,14,17,22]. 

destruction for both projectiles. Changing from C to Cu we increase the size 
of the targets, and therefore the probability for the halo neutrons being spec- 
tators decreases. Again the smaller size of 6 He is responsible for the larger 
difference between the long-dashed and short-dashed curves compared to the 
case of n Li projectile. The same behavior is observed for a Pb target. As the 
target mass increases the probability for the two neutrons being spectators 
decreases. The dominant process for a heavy target is the one where all the 
three halo constituents interact with the target. 



In Fig. 6 we plot interaction cross sections for the two projectiles and the 
three targets that we are considering. They are cross sections corresponding 
to processes where the projectile undergoes some reaction, and they are given 
by the sum of the two-neutron removal and core breakup cross sections. As 
a general rule <t_ c is larger than <7_2n, and only for low beam energies and 
heavy targets <7_ 2n dominates due to the large distance Coulomb collisions 
maintaining the core intact. Again, for 11 Li on Pb the computed interaction 
and two-neutron removal cross sections differ somewhat when the adiabatic 
cutoff energy huo a is increased from 0.15 MeV to 0.30 MeV. 

The agreement with the experimental data is good, except for the low 
beam energy point for the n Li projectile. However, this experimental point 
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Fig. 7. Interaction and two-neutron removal cross sections for 6 He and n Li on C, 
Cu and Pb for different values of the cutoff parameters. 

should be taken with caution, since the a-2n values for C and Pb targets 
are not consistent with the data for 11 Li on beryllium or gold, for which the 
experimental values of <r_ 2 n at 30 MeV/nucleon are 0.47±0.10 b and 5.0 ±0.8 
b, respectively [62]. 

As shown in Table 3 the values of the cutoff parameters R^j t are not 
uniquely determined. In Fig. 7 we examine the dependence of the interaction 
and two-neutron removal cross sections on these parameters. For a carbon 
target the main uncertainty appears in the case of a neutron spectator, where 
the value of R^ t can vary between 2.6 fm and 3.9 fm. We used 3.5 fm in the 
previous calculations. In Fig. 7 we compare the cross sections with R^ t —2.5 
fm and 3.5 fm, while we maintain 4.0 fm and 4.8 fm for R^Jt^ and R^\ 
respectively. The difference between these two calculations is very small. For 
Cu and Pb targets the largest range of variation for the cutoff parameters 
is for the core spectator. In the calculations shown up to now we have cho- 
sen the value obtained with procedure c), but procedures a) and b) can give 
significantly larger values for these parameters, see Table 3. 

In Fig. 7 we also compare aj and cr_2 n obtained with R\.^ = 5.7 fm and 
7.0 fm and Ri^ = 6.1 fm and 8.0 fm for a Cu target, and R^^ = 7.6 fm and 
10.0 fm and R cut = 7.9 fm and 10.0 fm for a Pb target. It is then striking that 
even when the difference between the parameters is substantial, the difference 
between the computed cross sections is rather insignificant. The main reason 
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Fig. 8. Longitudinal (solid) and transverse (long-dashed) core momentum distribu- 
tions after fragmentation of 6 He and n Li on C, Cu and Pb. The short-dashed line 
shows the contribution from core participant. Experimental data are from [9,16] for 
6 He on C, from [12,63] for n Li on C and from [12] for n Li on Pb. The experimental 
data have been scaled to the calculations. 

is that for these two targets an important contribution comes from the large 
impact parameter part of the Coulomb contribution (the reaction in Fig. 3d), 
where the cutoff parameters are unimportant. This is especially true for the 
lead target with the large Coulomb interaction. 

The main conclusion obtained from Fig. 7 is that the value chosen for 
the cutoff parameters within the range shown in Table 3 is not substantially 
modifying the cross sections. We may view these parameters as relevant for 
fine tune, while the cross section sizes and general behavior are maintained. 



4-4 Momentum distributions 



In Fig. 8 we show the longitudinal and transverse core momentum distri- 
butions after fragmentation of 6 He and n Li on on C, Cu and Pb. The calcula- 
tions have been done at a beam energy of 300 MeV/nucleon. For this energy 
the difference between longitudinal and transverse momentum distributions 
is small. The transverse distributions are in general wider than the longitudi- 
nal, typically by 6 to 12 MeV/c. The experimental data are taken from [9,16] 
for 6 He on C, and they are obtained at a beam energy of 400 MeV/nucleon 
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and 240 MeV/nucleon, respectively. For n Li on C the experimental data are 
taken from [12,63] and they have been obtained with a beam energy of 280 
MeV/nucleon. For n Li on Pb the data are from [12] and they correspond to a 
beam energy of 280 MeV/nucleon. Within this range of energies the difference 
in shape between the computed distributions is not visible. In all the cases the 
experimental data are given in arbitrary units, and they have been scaled to 
the calculations. The agreement between them and the calculations is good. 

In Fig. 8 the short-dashed line shows the contribution from core partici- 
pant. For light targets the dominating processes are those where the core is 
not participant, and therefore its final momentum is similar to the one that 
it had inside the halo projectile (except for the final neutron-core interaction, 
that does not modify significantly the core momentum [35]). When the target 
size increases the weight of the core participant contribution becomes more 
important. For Pb this contribution is very close to the total. The dominat- 
ing process is then the one shown in Fig. 3d, where only the low momentum 
transfer part of the Coulomb interaction should be included. The momentum 
transferred to the core in the collision is small compared to its initial momen- 
tum distribution, which therefore for a entirely different reason again is left 
relatively unchanged. Although in this case the three-body continuum wave 
function should be used, only the final neutron-neutron interaction has been 
included. Nevertheless, due to its large mass, core momentum distributions 
are not expected to be strongly modified by the final state interaction [34,35]. 
The consequence is that the shape of the core momentum distributions is al- 
most independent of the target. Only the initial core momentum distribution 
within the projectile is decisive. 

In Fig. 9 we show longitudinal and transverse neutron momentum dis- 
tributions after two-neutron removal fragmentation of 6 He and n Li on C, 
Cu and Pb. The calculations have been performed at a beam energy of 300 
MeV/nucleon. The difference between longitudinal and transverse distribu- 
tions is very small. Usually the computed transverse distributions are from 2 
to 4 MeV/c wider, and the difference between them is hardly visible. 

As in Fig. 8 the short-dashed curve shows the contribution from core- 
participant processes. The behavior is similar to the one found in Fig. 8. For 
light targets the core spectator contribution is dominating, while for heavy 
targets the main contribution comes from processes where the core partici- 
pates. The only available experimental neutron momentum distribution is for 
6 He on C [16] obtained at a beam energy of 240 MeV/nucleon. The agree- 
ment with the experiment is excellent, and only in the tails the computed 
curve is above the experimental distribution. This is due to the experimental 
neutron acceptance, that is limited in horizontal and vertical directions to a 
momentum of 50 MeV/c. 
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Fig. 9. The same as in Fig. 8 for neutron momentum distributions. The experimental 
data for 6 He on C are from [16]. 

For a lead target the neutron momentum distribution is dominated by 
core participant processes, and in particular by processes where the momen- 
tum transferred to the core is small. As discussed in Section 2, the final state 
is properly described for relatively small impact parameters by the reaction 
scheme shown in Fig. 1. For heavy targets the reaction is dominated by pro- 
cesses where mainly the low momentum transfer part of the Coulomb in- 
teraction contributes. The breakup process is then relatively gentle, and the 
three-body continuum wave function should be used to describe the final state 
as shown in Fig. 2. Therefore Eqs.(8) and (18) should then be used to describe 
these processes. 

It is well known that neutron momentum distributions are highly influ- 
enced by the final state interaction between fragments, while this influence 
is less important for core momentum distributions [34,35] due to the larger 
mass. In Figs. 8 and 9 we computed all distributions using the reaction scheme 
of Fig. 1, and therefore for the case of Pb target, where the final three-body 
continuum wave function should be used, only the final neutron-neutron in- 
teraction is included. As mentioned this has only a small influence on the core 
momentum distributions, where the agreement with the experiment shown in 
Fig. 8 is satisfactory. However, the computed neutron momentum distributions 
after fragmentation on a heavy target like Pb shown in Fig. 9 are expected to 
be too broad compared to the experimental data due to the neglect of the final 
state interaction between the neutron and the core for the dominating process. 
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Fig. 10. Various distributions obtained after fragmentation of He and Li on C 
at 300 MeV/nucleon. We only include the dominating parts where the a particle 
and one neutron are spectators, (a) is the transverse distribution of the 5 He center 
of mass momentum, (b) is the invariant mass spectrum of 5 He, (c) is the angular 
distribution with 6 as the angle between the relative a-neutron momentum and the 
center of mass momentum of 5 He, and (d) is the the same as (c) for a 11 Li projectile. 
The computed distributions have been convoluted with the instrumental response 
[65]. The experimental data are from [16] in (a) and (b), from [15] in (c) and from 
[66] in (d). 

These computed results, including the neutron-neutron final state interaction, 
are still much more realistic than the Fourier transform of the initial wave 
function implying that all final state interactions are neglected. 

For completeness we show in Fig. 10 four other interesting observables ob- 
tained after fragmentation of 6 He and 11 Li. As discussed above, inclusion of 
all the final state interactions is required for heavy targets corresponding to 
reaction scheme Fig. 2. Thus we focus here only on a C target. Similar ob- 
servables for the n Li projectile are already published in [20]. The calculations 
are performed for a beam energy of 300 MeV/nucleon, and the experimental 
data in Figs. 10a and 10b from [16] and in Fig. 10c from [15] are given in arbi- 
trary units for a beam energy of 240 MeV/nucleon. The ones in Fig. lOd from 
[66] are also given in arbitrary units for a beam energy of 287 MeV/nucleon. 
In ref. [64], where angular distributions for 6 He projectile are shown, we in- 
advertently used a wrong relative phase between interfering two-body final 
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states. The conclusions in [64] still hold but the numerical results now com- 
pare slightly better with measured values. 

The calculations in Fig. 10 are in agreement with the experiments, except 
for the angular distribution in Fig. 10c. The main reason for this discrepancy is 
due to the approximation of spheres instead of the proper cylinder geometries 
shown in Figs. 1 and 2. The exceptional sensitivity arises in this case, because 
the dependence on cos of the components with neutron-core relative p-wave 
angular momentum projection m = ±1 are very flat while the m = com- 
ponent is strongly varying [64]. Unlike the other quantities the proper cutoff 
symmetry is then decisive for this particular observable. The cylinder would 
remove relatively more of the m = component of the neutron-core p-wave 
and the computation would approach the measured distribution. The angular 
distributions are very different for 6 He and n Li projectiles as seen by compar- 
ing Figs. 10c and lOd. The large average value for 10 Li reflects the dominating 
constant contribution from neutron-core relative s-waves and the asymmetric 
variation is due to the interference of s and p-waves. The strong dependence 
on the proper cutoff of the different angular momentum projection states is 
then not present. It is especially interesting that the angular distribution for 
10 Li in Fig. lOd agrees so well since these computations were performed before 
the measurements were available. 



5 Summary and conclusions 

In this paper we formulate a model to describe fragmentation reactions 
of two-neutron halo nuclei on any nuclear target, from light to heavy. It is 
well established that nuclear and Coulomb interactions dominate for light and 
heavy targets, respectively. Thus both must be incorporated in a compre- 
hensive description which also should include intermediate target masses. We 
confine ourselves to relatively high energies where the reaction time can be 
expected to be short compared to the time scale of the intrinsic halo motion. 

The structure computations proceed in two steps, i.e. few-body models 
are used to describe the halo degrees of freedom and then the effects of the 
intrinsic structure of the halo constituents are considered. We adopt the same 
two-step strategy and start with a reaction model for weakly bound inert halo 
particles. Due to the spatial extension of the initial three-body system the 
dominating processes for short-range interactions are independent reactions 
of each halo particle (participant) with the target as for free particles. The 
participant is then elastically scattered or "something else" called absorption 
happens while the other two halo particles (spectators) continue undisturbed 
under the influence of their mutual interaction. 
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We use the optical model to describe the participant-target interaction. 
Elastic scattering is then described in details while different processes are 
considered as absorption which means that all other channels only remove 
probability and otherwise forgotten. However, this is the essence of the infor- 
mation in breakup experiments where the forward fragmentation products are 
measured. Furthermore, absorption in this optical model sense accounts for 
intrinsic degrees of freedom related to the target and the halo constituents. 
This is in other words consistent with the use of a three-body model for the 
halo structure. Going beyond this approximation in the reaction description 
requires at least as accurate a treatment of core degrees of freedom of the 
initial state. 

Another crucial ingredient to understand the reaction mechanisms is the 
finite sizes of both halo constituents and target. Conceptually this is related 
to the necessity both of accounting for intrinsic degrees of freedom and for de- 
scribing more than one simultaneous halo particle interaction with the target. 
This becomes clear in the impact parameter picture where more than one par- 
ticle at the same time may arrive at the target position within the interaction 
cylinder. We characterize the different processes of absorption and/or elastic 
scattering of the different halo particles that simultaneously interact with the 
target. For one of them (the core if possible) we employ optical potentials, 
and for the other(s) we use the part of the three-body wave function where 
this constituent is close enough to the first one such that it can be considered 
to be inside its cylinder. The model provides absolute cross sections for all 
different processes. 

So far we discussed the short-range parts of the interactions. In the im- 
pact parameter picture we may more conveniently formulate this as a re- 
striction to relatively large momentum transfer to the participant. Then the 
reaction mechanism is basically independent halo particle reactions dictating 
the division into two independent two-body final state systems (spectators 
and participant-target) and in our formulation resulting in the participant- 
spectator model. This division is not for convenience of computation but to 
include the proper physics. 

The large impact parameters or equivalently the small momentum trans- 
fers are important for long-range interactions. For nuclear halos this means the 
Coulomb interaction which is active when the charged core is a participant. 
When the target is heavy with a high charge this gives rise to the main contri- 
bution to the breakup cross sections arising from very distant collisions. The 
momentum transfer to the core then excites low lying states of the halo which 
for Borromean systems imply the low lying continuum states. These collisions 
are gentle and none of the halo particles are absorbed. Consequently the final 
state must consist of the three-body halo system and an independent target. 
The model distinguishes for these reasons between small and large momentum 
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transfer and allows therefore simultaneous treatment of nuclear and Coulomb 
interactions including interference terms. 

We use the model to investigate fragmentation reactions of 6 He and n Li 
on different targets, light as carbon, intermediate as copper, and heavy as 
lead. Total interaction, two-neutron removal and core breakup cross sections 
are computed as function of beam energy and found to be in good agree- 
ment with the available experimental data. We compute not only the terms 
arising from the dominating reaction mechanism, but also the parts usually 
neglected: Coulomb contribution for light targets and nuclear contribution for 
heavy targets. Reaction cross sections with targets of intermediate mass and 
charge, where nuclear and Coulomb interaction interfere and give comparable 
contributions, are also calculated. The model parameters are obtained from in- 
dependent sources like nuclear radii, binding energies and two-body scattering 
experiments. 

The cross sections for reaching the different final states consisting of the 
non-absorbed particles vary with target, projectile and beam energy. For light 
targets two neutrons and no core is most probable while all three halo particles 
or the core alone is least probable. For heavy targets absorption of all three 
particles is most probable except at low beam energies where simultaneous 
survival of all three halo particles is most probable. The smallest probability 
is found for core survival alone or the core with one neutron. 

Core destruction is more probable than two-neutron removal for all ener- 
gies and targets. Core destruction receives the largest contribution for heavy 
targets when all halo particles are absorbed while for light targets the two- 
neutron spectator contribution is comparable in size. Two-neutron removal re- 
ceives for heavy targets the largest contributions from core participation while 
for light targets the largest contribution arises from the core as a spectator. 
The relative sizes of all these cross sections reflect the reaction mechanisms 
which in this way are open for experimental tests. 

We have also computed core and neutron momentum distributions after 
fragmentation. This is a strong test of validity of the model, since total values 
obtained after integration could be in good agreement with experiments even 
when the non-integrated differential cross sections differ from measurements. 

It is known that final state interactions play an essential role in the momen- 
tum distributions, especially for light particles as neutrons. For light targets 
the dominant reactions are the ones with one or two spectators, and the model 
accounts properly for the final state interaction between spectators. The neu- 
tron momentum distribution receives a narrow contribution from one-neutron 
absorption and a broader contribution from neutron participation and sub- 
sequent scattering. The agreement with the available experimental core and 
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neutron momentum distributions is excellent for both projectiles. 

For heavy targets the low momentum transfer reaction dominates, and the 
appropriate reaction mechanism involves the low lying three-body continuum 
excitations in the final state. This implies that all three initial three-body 
interactions are necessary in the final state. In the numerical examples we 
only included two at a time. For core momentum distributions, where the fi- 
nal state interactions are expected to be less relevant, we still obtain good 
agreement with the experiment. For neutron momentum distributions we in- 
clude the neutron-neutron interaction but not the neutron-core interaction. 
The distributions are therefore expected to be somewhat too broad. 

The model predicts comparable shapes for core and neutron momentum 
distributions for light and heavy targets although the reaction mechanisms 
are very different. For light targets one- neutron absorption is the dominating 
process leaving the other neutron and the core untouched in the final state 
under their mutual interaction. For heavy targets the dominating process is a 
distant collision of low momentum transfer to the core leaving the two neutrons 
untouched and the core with a momentum very close to its initial value. Still 
the three particles are influenced by the final state interactions. In both cases 
the core and the neutrons would appear basically with their initial momentum 
distributions modified by final state interactions. This is not tested in details, 
but still in qualitative agreement with experiments. 

The model describes a number of different reaction mechanisms and pro- 
vides the appropriate relative weights for the corresponding processes. The 
same consistent model with one set of parameters is used throughout for all 
observables. In this paper the model is tested in details by comprehensive cal- 
culations for the most studied two-neutron halo nuclei. This includes absolute 
values of the differential cross sections as functions of beam energy and tar- 
get. Dominating contributions to many relative distributions can sometimes 
be explained by simpler models since they are surface or tail dominated. How- 
ever, the intrinsic degrees of freedom are necessary to obtain reliable absolute 
values, especially for less dominating processes. The details are unnecessary 
for breakup reactions as long as predictable total removal probabilities can be 
obtained as provided by the phenomenological optical model. 

In conclusion, the model seems to be very successful but the many pre- 
dictions should be tested by experiments to tie down in details the correct 
reaction mechanisms. 
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